A simple derivation of the Gompertz law for human mortality. 
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The Gompertz law for the human mortahty rate as a function of age is derived from a simple 
model of death as a result of an exponentially rare escape of abnormal cells from the immunological 
response. 



Human life is finite as are the lives of some well known 
physical objects for example radioactive nuclei. The 
probability, S{t), that a given nucleus will survive time t 
is 

S{t) = exp(-i/io). (1) 

If survival of humans were governed by the same law 
with, for example, to = 70 years we would have millions 
of people with the age, say, Stg — 350 years. However, 
the longest recorded human life lasted only 122 years. 
Apparently, the law of human mortality should be dras- 
tically different from Eq. 

The famous astronomer Halley^ and the great math- 
ematician EulerA were the first to attack this problem, 
but the law was empirically found only in 1825 by the 
actuary and self-taught mathematician Benjamin Gom- 
pert»3. If S(t) is probability of a human surviving till 
age t years then the mortality function is defined as 

fi{t) = -d\nS{t)/dt. (2) 

Gompertz found that when childhood diseases are over- 
come {t > 25) the statistics of mortality obeys the law 

/i(t) =A*(25)exp([i-25]Ai), (3) 

where ti ~ 10. This law should be contrasted with the 
constant mortality function fi{t) = I/Iq following from 
Eq. The Gompertz law means that S{t) decays dou- 
ble exponentially, much faster than Eq. practically 
eliminating people older than 122 in the world popula- 
tion. Statistical data for many countries and three cen- 
turies confirm the Gompertz law in the range of 
covering more than three orders of magnitude (see Refs.** 
and references therein). The law works also for different 
species such as rats, mice and fruit files. Thus, the Gom- 
pertz law emerges as one of the greatest quantitative laws 
of biology. 

There were a number of attempts to derive the Gom- 
pertz law. The most quoted derivation uses the lan- 
gauge of the reliability theory designed for man-made 
machines^ and is quite complicated. Here, we would 
like to suggest a simple version of derivation of Eq. 
based on a naive understanding of immunology. Think- 
ing about cells, we assume that a population of defective 
cells (mutated cells that do not fulfill their normal func- 
tion, cancer cells undergoing uncontrolled proliferation, 
cells that produce masses of defective protein forming 



amyloids and senile plaques in aged tissues leading to 
diseases like Alzheimer) becomes fatal if the normal or- 
ganism defense system gives this population a time r 
necessary to reach a critical size. Let us assume that at 
age 25, bad cells would encounter strong and fatal im- 
mune response during time r on average A'o 3> 1 times. 
(We can imagine a bad cell as criminal who is being elim- 
inated at any encounter with randomly patrolling police- 
men (macrophages, natural killers, or apoptosis), before 
growing strong enough to successfully resist them. Until 
this time, on average, he would meet iVo policemen.) For 
absolutely random encounters, the probability P{N) of 
the number N of encounters with immune response sys- 
tem during time t is given by the Poisson distribution, 

P(N) ^ ^^f^Ej^, (4) 

where A^o is the average value of N. The probabil- 
ity of a population of defective cells reaching the crit- 
ical size and causing the death of the host at age 25 
is ^(25) = P(0) = exp(-Ao). It is known that ^(25) 
is of the order of 3 10"^. This means that, indeed, 
Nq ~ 8 ^ 1. Let us assume now that the immune re- 
sponse slowly weakens with the age so that at t > 25 the 
average number of encounters during " microscopic" time 
T decreases linearly as No{t) = No — {t — 25)/ti. This may 
happen due to accumulation of mutations in immune re- 
sponse cells, their limited potential for self-renewal, or 
an overall decay of the organism energetics. (In the lan- 
guage of criminals and policemen, this would mean slow 
decay of the number of patrolling policemen, for exam- 
ple, due to budget restrictions). Substituting this No{t) 
instead of Nq into the Poisson formula for P(0) we arrive 
at fi{t) — exp[— Ao(t)] and, therefore, at Eq. 

Notice that a relatively small change of No{t) leads 
to a very strong exponential growth of /i(t). A linear 
decay of A'o (t) can be the initial stage of an exponential 
relaxation No{t) = Noexp{—[t — 25]/Noti). In this case, 
at large ages t ~ 100 we arrive at downward deviations 
of /i(t) from the Gompertz law. This qualitatively agrees 
with somewhat slower than Eq. (jSJ growth of mortality at 
t > 100 cited indi^i^. These empirical deviations from 
Eq. ^ are still under discussion because of relatively 
poor statistics for t > 100. 

I am grateful to M. Azbel, A. Finkelstein, A. 
Chklovski, D. Chklovskii, S. Grigoryev, I. Ruzin, L. 
Shklovskii and D. Stauffer for useful discussions. 



2 



^ E. Halley, Phil. Trans. Roy. Soc. 17, 596 (1693). 

^ L. Euler, Histoire de I'Academie Royale des Sciences et Belle 

- Letters, 144 (1760). 
^ B. Gompertz, Philos. Trans. Roy. Soc. Lond. A 115, 513 

(1825). 

* D. Stauffer, The Complexity of Biological Ageing, p. 131 in: 
Thinking in Patterns, ed. by M.M. Novak, World Scientific, 



Singapore 2004; |cond-mat/0310038l 

L. A. Gavrilov, N. S Gavrilova, J. Theor. Biology 213, 527 
(2001). 

^ M. Azbel, Physica A 249, 472 (1998). 
R. M. C. de Almeida, S. Moss de Oliveira, and T. J. P. 
Penna, Physica A 253, 366 (1998). 



